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ABSTRACT 

An  experimental  study  was  performed  to  determine  the  feasibility  of 
a  laboratory  experiment  to  test  an  existing  theoretical  model  describing 
sound  propagation  into  a  fast  bottom  underlying  a  wedge  shaped  medium. 
Sand  under  fresh  water  was  found  to  satisfy  the  constraints  of  the  theor- 
etical model  and  to  simulate  the  continental  shelf.  In  a  laboratory  ex- 
periment, accuracy  of  density,  sound  speeds,  and  attenuation  was  shown 
to  be  sufficient  to  allow  quantitative  comparison  to  the  predictions  of 
the  beam  angle. 
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I.   INTRODUCTION 

The  behavior  of  radiation  in  a  wedge-shaped  medium  overlying  a  fast- 
bottom  medium  has  been  investigated  both  theoretically  and  experimentally 
with  both  optical  and  acoustic  sources.  For  example,  in  1971  Tien  and 
Martin  [Ref.  1]  examined  the  behavior  of  a  laser  beam  coupled  into  a  thin, 
tapered,  dielectric  film  deposited  on  a  substrate  with  a  higher  refrac- 
tive index.  As  the  light  propagated  toward  the  apex  of  the  wedge  shaped 
film,  perfect  reflection  was  observed  until  the  changing  angle  of  inci- 
dence decreased  below  the  critical  angle,  then  the  light  was  converted 
into  radiation  in  the  substrate. 

An  acoustic  analog  to  the  above  optical  study  was  performed  by 
Kuznetzov  [Ref.  2]  in  1973.  Kuznetzov  developed  a  normal  mode  theory 
for  sound  propagation  in  both  a  wedge  and  underlying  half-space  substrate, 
where  the  substrate  sound  speed  was  faster  than  that  of  the  wedge. 
Kuznetzov  concluded  that: 

1.  Sound  traveling  toward  the  wedge  apex  would  be  total  1  reflected 
until  the  angle  of  incidence  decreased  to  the  limiting  angle  of  total 
reflection. 

2.  Any  sound  incident  at  less  than  the  limiting  angle  of  total  re- 
flection would  be  completely  refracted  into  the  underlying  half-space. 
This  total  refraction  would  occur  along  the  wedge/half-SDace  boundary 
from  the  wedge  apex  to  the  point  where  the  limiting  angle  was  first 
achieved. 

3.  Acoustic  energy  in  the  half-space  would  be  columnated  into  a 
well-defined  beam.  The  beam's  maximum  pressure  would  occur  at  an  angle 
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of  depression  (measured  from  the  plane  of  the  wedge/half-space  interface) 
that  would  lie  between  8  and  28,  where  8  is  the  wedge  angle. 
Kuznetzov  performed  a  series  of  experiments  that  supported  his  theory. 

The  formation  of  a  well-defined  beam  within  the  substrate  (for  the 
optical  case)  was  subsequently  demonstrated  by  Tien,  Smolinsky  and 
Martin  [Ref.  3].  Reference  3  also  presented  two  theoretical  results: 
(1)  ray-optics  predicted  refraction  into  the  substrate  beyond  the  cutoff 
distance  (the  point  at  which  the  critical  angle  was  reached)  and  (2) 
theory  predicted  refraction  into  the  substrate  before  the  cutoff  distance, 

A  similar  well-defined  beam  was  predicted  and  observed  by  Sigelman, 
et  al .  [Ref.  4  J  in  a  water-aluminum  system.  Maximum  pressure  occurred 
at  11  degrees  below  the  water-aluminum  interface  when  the  wedge  angle  was 
1.3  degrees.  This  was  well  outside  the  range  predicted  by  Kuznetzov. 
The  observed  beam  was  broader  than  predicted. 

In  1980,  Bradshaw  [Ref.  5]  extended  a  computer  model  first  developed 
by  Kawamura  and  Ioannou  [Ref.  6j  to  describe  the  behavior  of  sound  in  a 
fast  bottom  underlying  a  wedge  shaped  medium.  The  model  predicted  the 
formation  of  well-defined  beams  in  the  bottom,  and  predicted  the  effect 
of  attenuation  in  the  bottom  on  the  beam.  A  comparison  of  the  predic- 
tions of  this  model  with  the  experimental  data  of  Netzorg  [Ref.  7]  gave 
qualitative  agreement,  but  since  the  experiment  did  not  match  all  model 
constraints,  quantitative  comparison  was  impossible.  It  is  the  purpose 
of  this  laboratory  analysis  to  determine  if  a  laboratory  experiment 
that  fulfills  the  constraints  of  the  model  can  be  designed.  If  at  all 
possible,  it  is  further  desired  to  model  the  continental  shelf  in  the 
design  of  the  experiment. 


II.  CONSTRAINTS 

The  general  design  of  the  proposed  experiment  is  illustrated  in 
Figure  1.  Constraints  will  be  imposed  by  both  the  theoretical  model 
and  the  measurements. 

A.  MODEL  CONSTRAINTS 

The  Bradshaw  model  [Ref.  5 J  requires  that  the  sound  speed  in  the 
bottom  exceed  that  in  the  overlying  wedge.  Physical  parameters  required 
to  be  known  are  the  density  and  sound  speed  of  the  two  media,  the  slope 
of  the  wedge,  and  the  attenuation  of  the  substrate  medium.  A  further 
constraint,  dependent  on  the  substrate  sound  speed,  will  be  discussed  in 
detail  later.  This  model  assumes  planar  waves  incident  on  the  wedge- 
bottom  interface. 

3.  OTHER  CONSTRAINTS 

Since  it  was  desired  to  model  a  real -world  environment,  specifically 
the  continental  shelf,  a  laboratory  set-up  consisting  of  fresh  water 
over  sand  was  an  obvious  first  choice.  The  sound  speed  of  fresh  water 
was  expected  to  be  about  1481  m/s  [Ref.  8]  and  the  speed  in  sand  about 
1700  m/s  [Ref.  9j.  Thus,  the  requirement  of  a  faster  bottom  medium 
would  be  readily  fulfilled. 

A  wedge  of  water  would  require  a  high  frequency  (i.e.,  greater  than 
100  kHz)  to  permit  planar  incident  waves  in  a  laboratory-sized  experi- 
ment. If  one  hundred  wavelengths  are  desired  to  ensure  plane  waves,  at 
100  kHz  the  source  would  need  to  be  1.5  m  away  from  the  wedge  apex. 
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The  use  of  a  water/sand  system  would  permit  the  hydrophone  to  be 
moved  about  in  both  media;  the  use  of  plaster  to  simulate  rock,  for 
example,  would  not.  While  it  would  also  be  necessary  to  maintain  a  flat 
slope  of  at  least  five  degrees  between  the  two  media,  a  water-sand  sys- 
tem was  expected  to  be  adequately  stable.  The  smoothness  of  the  slope 
is  required  by  the  model;  the  five  degree  slope  encompasses  the  maximum 
for  the  continental  shelf. 
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Figure  1.   Model  Geometry 
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III.  EXPERIMENTAL  DESIGN 

A.  MATERIAL  SELECTED 

Fresh  (tap)  water  and  #30  fine  sand  were  the  media  used  in  the  experiment 
The  grain  size  of  #30  fine  sand  varies  from  0.70  mm  to  0.15  mm.   It  was 
readily  demonstrated  that  the  sand  could  maintain  a  slope  as  large  as  25 
degrees.  To  remove  all  air  bubbles,  the  water/sand  system  was  heated 
with  immersion  heaters.  Bleach  was  added  at  a  ratio  of  one-half  gallon 
bleach  to  70  gallons  of  water  plus  sand  in  order  to  control  biologic 
growth. 


B.  VELOCITY  MEASUREMENT  TECHNIQUES 

An  excellent  summary  of  methods  for  measuring  sound  speeds  in  rocks 
and  minerals  is  given  by  Anderson  and  Liebermann  [Ref.  10J.  Specialized 
methods  for  use  on  small  specimens,  iron  compounds,  etc.  or  for  high 


nated  as  inappropriate  for 
c  interferometry  was  also  dis- 
01  to  0.04  percent)  as  it  is 
zes  frequencies  (10-60  MHz) 


pressure  and  temperature  systems  were  elimi 
the  experiment  being  considered.  Ultrasoni 
carded  despite  its  exceptional  accuracy  (0. 
best  applied  to  nonporous  samples  and  utili 
that  would  cause  unacceptable  attenuation  in  the  verification  experi- 
ments. Resonance  methods  with  appropriate  frequencies  and  accuracies 
were  undesirable  as  they  are  applicable  to  crystalline  specimens  only 
Ultimately,  pulse  transmissions  were  selected  as  they  best  fit  the 
necessary  criteria.  This  method  is  applicable  to  either  fluid  or  por 
ous  materials,  acceptable  over  a  frequency  range  of  50  kHz  to  10  MHz, 
and  accurate  to  within  one  to  three  percent. 
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1 .  Size  Criteria 

Steel -bound  glass  tanks  measuring  70  cm  square  by  60  cm  high 
were  readily  available.  Length  vs  depth  criteria,  from  Reference  10, 
require  that  the  layer  of  sand  on  the  tank  bottom  should  not  be  less 
than  14  cm  to  avoid  boundary  reflection  problems.  To  satisfy  this  con- 
straint, a  sand  layer  of  20  cm  was  chosen. 

To  avoid  pulse  dispersion,  Reference  10  also  suggested  that  the 
sand  depth  not  be  less  than  five  wavelengths.  Since  the  sound  speed  of 
the  sand  was  expected  to  be  1700  m/s,  and  the  frequency  would  be  of  the 
order  of  100  kHz,  five  wavelengths  would  be  8.5  cm.  A  sand  depth  of 
20  cm  would  also  satisfy  this  restraint. 

To  avoid  particulate  scattering,  the  wavelength  must  be  at 
least  three  times  the  grain  size  [Ref.  10].  Since  the  largest  grain 
size  in  number  30  sand  is  0.070  cm,  a  wavelength,  equal  to  1.7  cm, 
readily  qual ifies. 

Since  sound  speed  is  not  a  function  of  frequency  [Ref.  9],  re- 
sults from  this  laboratory  analysis  were  expected  to  be  applicable  to 
frequencies  used  in  the  experiments. 

2.  Electronic  Equipment 

A  schematic  of  the  equipment  configuration  is  shown  in  Figure  2. 
All  components  were  off-the-shelf  and  readily  available. 

Output  from  a  General  Radio  model  1310  oscillator  with  a  fre- 
quency range  of  2  Hz  to  2  MHz  was  fed  simultaneously  into  a  frequency 
counter  and  a  tone  burst  generator.  The  counter,  a  Hewlett-Packard 
5233L,  would  read  +10  Hz  at  100  kHz.  The  General  Radio  Type  396-A  tone 
burst  generator  was  used  to  generate  either  8  or  16  cycle  pulses.  The 
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output  passed  through  a  Hewlett-Packard  HP  467-A  power  amplifier  before  being 
fed  to  a  transducer.  The  pulsed  signal  was  also  fed  to  the  trigger  in- 
put of  a  Tektronix  type  Rm-503  oscilloscope. 

Three  types  of  transducers  were  used  as  sources.  The  first  was 
a  homemade,  mylar,  broad-band  transducer  with  a  3  cm  by  9  cm  active 
face.  This  source  was  quite  directional  at  100  kHz  but  proved  to  be 
difficult  to  move  through  the  sand  (for  sound  speed  and  attenuation 
measurements)  and  was  not  sufficiently  powerful  to  generate  usable  sig- 
nals in  the  sand  (as  would  be  required  in  model  verification  runs).  The 
mylar  transducer  required  by  150  V  DC  supply  and  polarizing  network. 

A  second,  smaller,  homemade  transducer  was  tried  as  a  source. 
This  transducer,  a  matrix  of  ceramic  elements,  has  a  1  x  2  cm  active 
face  and  was  highly  directional  at  100  kHz.  This  transducer  offered 
little  improvement  over  the  first. 

Celesco  Industries  type  LCI 0  hydrophones  had  been  consistently 
used  as  receivers.  The  LCI 0  is  a  small  (0.97  cm  diameter  by  2.87  cm 
long)  cylinder,  with  a  receiving  range  of  0.1  to  120,000  Hertz.  The 
transducer  is  designed  to  be  omnidirectional  in  a  plane  perpendicular 
to  the  axis  of  the  cylinder  with  a  tolerance  of  +1  dB  at  100  kHz,  and 
omnidirectional  in  a  plane  containing  the  axis  of  the  cylinder  with  a 
tolerance  of  +2  dB  at  25  Hz.  The  omnidirectional  characteristic  was 
not  desirable  for  the  source,  as  this  would  increase  the  number  of  pos- 
sible reflected  signals  from  the  walls,  water  surface,  etc.,  but  the 
small  size  and  higher  signal  strength  were  desirable.  Geometry  was  used 
to  isolate  desired  signals  and  will  be  described  when  appropriate. 
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Because  the  direct  and  reflected  pulses  arrived  from  different 
directions,  use  of  LClO's  as  the  source  and  receiver  for  attenuation 
measurements  initially  produced  non-reproducible  results.  Directivity 
measurements  in  the  plane  perpendicular  to  the  axis  of  the  cylinder  were 
made  and  the  results  are  shown  in  Figs.  3  and  4.  Subsequent  measurements 
used  LC10  number  2338  as  a  source  and  LC10  number  2319  as  a  receiver. 
The  receiver  was  oriented  to  use  the  area  about  50  degrees  as  the  receiv- 
ing face.  The  transducers  were  always  arranged  geometrically  to  use  the 
radial  plane. 

The  received  signals  were  amplified  20  dB  by  a  Hewlett-Packard 
HP-465A  amplifier,  then  passed  through  a  Spencer-Kennedy  Laboratories, 
Inc.  model  302  variable  electronic  filter  (set  at  60  kHz  high  pass)  to 
eliminate  low  frequency  mechanical  noise  present  in  the  laboratory  be- 
fore being  passed  to  the  oscilloscope.  The  data  of  Fig.  5  demonstrates 
that  spherical  spreading  (1/r)  was  observed  for  all  LC10  source-receiver 
separations  greater  than  8.4  cm.  All  measurements  were  made  under  far- 
field  conditions. 
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Figure  3.     LC10  2338  Directivity  Pattern 
(units  of  volts) 
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Figure  4.      LCIO  2319  Directivity  Pattern 
(units  of  volts) 
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IV.  MEASURED  DATA 

For  all  of  the  following  measurements,  the  sand  layer  in  the  tank 
was  kept  flat  and  horizontal,  as  only  physical  parameters  were  being  de- 
termined and  not  wedge  characteristics. 

A.  DENSITY 

1.  Water 

According  to  Lange's  "Handbook  of  Chemistry"  [Ref.  11],  the 

density  of  distilled  water  ranges  from  0.99913  g/cm3  at  15°C  to  0.99707 

3 
g/cm  at  25°C.  The  expected  density  of  room  temperature  water,  to  3 

3 
significant  figures,  was  therefore  1.00  g/cm  .  Measured  volumes  of  the 

bleached  water  were  periodically  weighed,  and  the  observed  density  was 

3 
1.00  g/cm  ,  with  no  noted  variation. 

2.  Sand 

The  density  of  water  saturated  sand  was  measured  by  partially 

filling  a  weighed  100  ml  graduated  cylinder  with  a  water  over  sand  mix- 

3 
ture,  observing  the  volume  of  each,  and  using  a  density  of  1.00  g/cm 

for  the  water  to  calculate  the  weight  of  the  water  volume.  The  density 

of  the  water  saturated  sand,  for  eight  separate  measurements,  was 

1.98+0.03  g/cm3. 

B.  SOUND  SPEED 
1.  Water 

Three  separate  sets  of  sound  speed  in  water  measurements  were 
made.  The  first  involved  the  use  of  the  mylar  transducer  as  the  source. 
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One  LCIO  was  placed  20.1  cm  away,  and  the  second  LC10  was  moved  along 
the  straight  line  defined  by  the  source  and  first  receiver.  The  time 
of  flight  between  the  two  receivers  was  measured.  A  simultaneous  check 
for  1/r  dependence  confirmed  these  measurements  were  all  in  the  far- 
field.  The  data  are  summarized  in  Table  1  of  Appendix  A.  The  average 
sound  speed  was  1454  m/s  with  an  (N-l)  variance  of  33  m/s. 

The  second  set  of  data  utilized  a  reflected  signal  to  increase 
the  total  path  length.  The  mylar  source  was  fixed,  and  the  receiver, 
placed  between  the  source  and  the  wall,  was  moved  along  a  perpendicular 
line  between  the  source  and  the  wall.  The  time  of  flight  was  the  time 
difference  between  the  arrival  of  the  direct  and  reflected  signals;  the 
distance  was  twice  the  distance  from  the  receiver  to  the  wall.  These  re- 
sults are  summarized  in  Table  2  of  Appendix  A.  The  average  speed  was 
1469  m/s,  with  a  variance  of  14.7  m/s. 

In  the  third  set  of  measurements,  the  distance  over  which  the 
time  of  flight  was  measured  was  fixed.  The  LC10  source  was  suspended 
11.0  cm  from  the  tank  wall,  the  second  LC10,  placed  further  from  the 
wall  than  the  source,  was  moved  along  the  perpendicular  formed  by  the 
wall  and  the  source.  As  the  receiver  was  moved,  the  distance  differ- 
ence between  the  direct  path  signal  and  the  signal  reflected  from  the 
wall  remained  a  constant  22.0  cm.  The  time  difference  between  the  ar- 
rival of  the  two  signals  was  recorded.  Table  3  of  Appendix  A  summarizes 
the  results.  The  average  speed  was  1460  m/s,  the  variance  was  19  m/s. 

The  data  compared  well  with  theory.  From  Kinsler  and  Frey 

[Ref.  8],  the  equation  of  sound  for  distilled  water  is  given  by 

c  =  1403  +  5t  -  0.06t2  +  0.0003t3  (1) 

w 
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where  t  is  the  temperature  in  degrees  Celcius.  A  comparison  of  measured 
and  theoretical  values  is  presented  in  Table  VI-1.  Note  that  direct 
path  results  showed  greater  variation,  but  all  results  were  within  two 
percent  of  expected  values. 
2.  Sand 

Five  separate  sets  of  sound  speed  measurements  were  made  in  the 
sand;  four  in  conjunction  with  attenuation  data,  and  one  to  check  for  a 
sound  speed  gradient.  All  three  source  transducers  were  used.  For  the 
attenuation  runs,  the  source  was  buried  at  a  fixed  depth,  probing  was 
with  an  LC10  receiver  inside  a  glass  tube  for  the  depth  of  maximum  signal 
amplitude.  The  glass  tubes  were  necessary  to  prevent  abrasive  damage  to 
the  receiver.  The  data  are  summarized  in  Tables  1  through  4  of  Appendix 
B. 

The  fifth  set  of  data  utilized  the  two  LClO's  fixed  30  cm  apart. 
The  transducers  were  initially  placed  in  water  and  sand  was  added  to 
progressively  bury  them.  The  entire  process  was  carried  out  underwater, 
so  degassed,  saturated  sand  was  used.  The  data  are  summarized  in  Table 
5  of  Appendix  B.  Table  IV-2  lists  the  resultant  sound  speeds. 

a.  Theoretical  Sound  Speed  in  Sand 

Since  a  search  of  the  literature  revealed  many  theoretical 
and  experimental  values  for  the  speed  of  sound  in  sand  under  brine,  but 
not  under  fresh  water,  and  since  the  values  given  were  greater  than 
those  observed  in  this  experiment,  the  following  analysis  was  performed 
to  derive  a  theoretical  value  for  sound  speed  in  sand  under  fresh  water. 

From  Urick  [Ref.  12],  the  density  of  a  liquid-sediment 
mixture  is  given  by 
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pmix  =  YPw  +  (1  "  Y)ps  (2) 

where  pmix  =  the  density  of  the  mixture, 
Y    =  the  porosity, 
Pw   =  density  of  water, 

and    Ps   =  the  density  of  individual  sand  grains. 

From  equation  (2),  using  the  measured  value  p  .  =1.98  g/cm3  and 

mix       r 
3 
Ps  =  2.69  g/cm  from  Reference  9,  gives  y  =   0.42.  Also  from  Reference 

12,  the  speed  of  sound  in  the  saturated  sediment  is  given  by 

cmix=  [<V  +  ps(1  -T))(kwY+  ks(l  -  Y))]"1/2     (3) 
where  kw  and  k$  are  the  compressibility  of  the  water  and  sand,  respec- 
tively. Values  for  k  and  k  were  calculated  from  the  relationships 

c„  =  (o  k  )"1/2  (4) 

w   XKw  w7  K    ' 

and  cs  =  (Psks)"1/2  (5) 

Setting  cw  =  1481  m/s  [ftef.  8],  and  p  =  1.0  g/cm  ,  and  solving  for 


kw  yields 


k  ,  =  4.56  x  10"13  m  sec2  g"1  (6) 

w 


The  sand  in  the  experiment  was  composed  of  quartz  particles.  From 
Anderson  and  Lieberman  [Ref.  10],  page  360,  c  was  set  to  7000  m/s. 
Solving  equation  (5)  yields 

k$  =  7.59  x  10'15  m  sec2  g"1  (7) 

Substituting  the  above  values  in  equation  (3)  and  solving  for  c  . 

yields  an  expected  speed  for  sand  under  fresh  water  of  1605  m/s. 

A  comparison  of  the  expected  value  with  the  sound  speeds 
summarized  in  Table  IV-2  shows  close  correlation  with  the  value  obtained 
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TABLE  IV-2.  Measured  Sound  Speeds  in  Sand 


Measurement 
Method 

Variable  distance 
Variable  distance 
Variable  distance 
Variable  distance 
Fixed  source/receiver 


Sound  Speed 
(m/s) 

(N-l)  a 
(m/s) 

1579 

58.6 

1597 

44.7 

1434 

71.9 

1567 

32.6 

1613 

0.0 
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using  a  fixed  source-receiver  distance  (1613  m/s)  but  not  close  agreement 
with  the  majority  of  values  obtained  from  measurements  utilizing  glass 
tubes.  The  use  of  glass  tubes  disturbs  the  medium,  and  the  distance  from 
source  to  receiver  is  not  accurately  measurable.  Over  the  short  distances 
used,  a  0.5  cm  error  in  the  distance  would  lead  to  significant  errors. 
Thus,  the  value  obtained  from  the  fixed  source  and  receiver  method  is  more 
likely  accurate,  and,  not  unexpectedly,  shows  closer  agreement  with  the 
theoretical  values. 

b.  Sound  Speed  Gradient 

As  illustrated  in  Fig.  6,  Hamilton  [Ref.  13]  predicts  a  pro- 
nounced sound  speed  gradient  in  a  fine  sand  sediment.  A  change  in  slope 
of  27  s"  is  specified  at  zero  depth.  Although  the  gradients  were  estab- 
lished for  sand  under  brine,  they  are  attributed  to  a  decrease  in  sediment 
porosity  with  pressure,  to  temperature  effects,  and  to  pressure  effects  on 
the  pore  water,  and  would  therefore  be  expected  to  be  present  in  fresh 
water  systems  as  well.  Sound  speed  measurements  tabulated  in  Table  5  of 
Appendix  B  were  taken  specifically  to  check  for  the  presence  of  a  gradient. 
None  was  found  within  the  accuracy  limitations  of  the  method.  For  a  sound 
speed  of  1613  m/s,  a  2%   error  (see  Sec.  Ill -A)  would  be  32  m/s  and  would 
mask  a  27  m/s  change.  Even  if  a  27  m/s  gradient  were  present,  the  re- 
sulting 62  m  radius  of  curvature  would  not  be  sufficient  to  be  detectable 
over  the  experimental  ranges.  The  lack  of  a  detectable  gradient,  although 
not  mirroring  actual  continental  shelf  conditions,  is  consistent  with  the 
isovelocity  substrate  assumed  in  Bradshaw's  model  [Ref.  5]. 
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reference  13.) 
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C.  ATTENUATION 

1 .  Model  Limitations 

It  was  mentioned  in  Section  II-A  that  a  further  constraint,  to  be 
defined  once  the  substrate  sound  speed  had  been  determined,  must  be  ful- 
filled for  verification  of  Bradshaw's  model.  This  constraint  is  that  the 
signal  must  propagate  beyond  a  specific  range  determined  as  follows.  The 
model  uses  a  Green's  function  analysis  to  obtain  the  pressure  and  phase 
distribution  of  the  sound  in  the  bottom.  To  solve  the  appropriate  Green's 
function,  it  was  necessary  to  assume  that 

k2|r  -  r*|  »  2tt  (8) 

where  k^  =  co/c  is  the  wave  number  for  the  sound  in  the  sand.  Using  the 
determined  sound  speed  of  1613  m/s  and  assuming  a  frequency  of  100  kHz, 
Y.~   equals  389  m"  .  Thus,  it  is  required  that  |r  -  r"|  »  0.016  m 

The  quantity  \r   -  V     is  defined  in  Fig.  7.  The  quantity  r'  is 
the  distance  from  the  wedge  apex  to  a  point  on  the  wedge-substrate 
boundary  (integration  in  the  model  is  carried  out  over  r").     The  quantity 
r  is  the  distance  from  a  point  along  the  wedge-substrate  boundary,  called 
X,  to  the  field  point.  The  distance  X  is  the  distance  at  which  the  cri- 
tical angle  is  reached  and  where  lowest  mode  transmission  into  the  sub- 
strate begins.  The  distance  is  defined  by 

A   4  f  sinec  sinB'  K    ' 

6r  =  arccos  (c/c.)  (10) 

where  f  is  the  frequency, 

6  is  the  wedge  angle 

c  ,  c  are  the  sound  speeds  in  the  wedge  fluid  and  substrate, 
w   s 

respectively, 
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and    9c  is  the  critical  angle. 

For  f  =  100  kHz,  a  wedge  angle  of  3°,  c  =  1469  m/s  and  c  =  1613  m/s,  X 

w  s 

would  equal  0.17  m. 

Typically,  within  the  model,  the  integration  limit  is  from  nX  to 
NX,  where  n  =  0,1,2,...,  N  =  1,2,3,...,  N>n.   In  the  simplest  case,  the 
integration  would  be  carried  out  from  0  to  X.  This  places  a  constraint 
on  r,  and  thus  on  permitted  attenuation. 

Consider,  as  a  minimum,  a  factor  of  10  to  satisfy  the  "much 
greater  than"  inequality;  this  would  require  \r  -  r"\    to  be  greater  than 
0.16  m.  This  implies  mathematically  that  no  measurements  may  be  made  at 
point  r  within  a  circle  of  radius  16  cm  about  the  point  r" ,   and  since  r' 
can  vary  from  0.0  to  0.17  m,  an  area  of  "forbidden  territory"  is  traced 
about  the  wedge  apex  and  boundary  surface.  Thus,  to  test  the  model, 
since  measurements  must  be  made  outside  this  area,  the  attenuation  must 
not  be  so  great  that  insufficient  signal  remains  beyond  the  forbidden  zone 
(If  "much  greater  than"  were  taken  to  be  100  times,  measurements  would 
have  to  be  made  over  1.6  m  from  the  boundary). 
2.  Results 

The  methodology  for  collecting  attenuation  data  was  described  in 
Section  IV-B-2.  The  data  were  processed  graphically.  Taking  the 
natural  logarithm  of  the  well-known  equation 

V  -ar 

V  =  -°e  (11) 

r 

where  V  is  the  measured  voltage, 

V  is  the  source  voltage, 

r  is  the  distance  from  the  source, 
and    a  is  the  attenuation, 
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the  linear  relation 


ln(Vr)  =  ln(VQ)  -  ar  (12) 


is  obtained.  Thus  a,  in  nepers  per  meter,  will  be  the  slope  of  a  graph 
of  ln(Vr)  vs  r.  Graphs  of  the  four  data  sets  of  Appendix  B  are  shown  in 
Figures  8  through  11.  A  summary  of  the  attenuation  values  is  presented 
in  Table  IV-3. 

Figure  12  shows  a  plot  of  the  attenuation  vs  frequency.  The  line 
indicates  the  slope  of  any  line  having  a  dependence  of  attenuation  on  the 
first  power  of  frequency.  All  data  fall  within  the  limits  of  data  col- 
lected by  Hamilton  [Ref.  14J. 

3.  Reflection  Coefficients  and  Signal  Level 

Normal  incidence  pressure  reflection  and  transmission  coefficients 
were  measured  as  an  accuracy  check  on  the  measured  p  and  c  values,  and 
also  to  provide  an  estimate  of  the  sound  pressure  level  that  would  be 
transmitted  from  the  wedge  to  the  underlying  medium.  A  summary  of  the 
measured  data  is  presented  in  Appendix  C. 

The  reflection  coefficient  R  is  calculated  from 

R  =  !pwcw-  pscs>  (13) 

lpwcw+  pscsl 

and  the  transmission  coefficient  T  from 

T  =  1  -  R  (14) 

With  p  =1.00  g/cm3,  p„  =  1.98  g/cm3,  c  =  1469  m/s,  and  c  =  1613  m/s, 

vv  s  "  3 

R  is  calculated  to  be  0.37  and  T  to  be  0.63.  The  measured  values  of  re- 
flection were  0.31  and  0.29.  The  measured  transmission  coefficient  was 
0.65.  These  values  indicate  that  the  measured  density  and  sound  speed 
values  are  within  the  expected  range. 
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Figure  12 .  Attenuation  vs  Frequency 
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From  Table  1  of  Appendix  C,  a  received  voltage  of  6  volts  is 
typical.  Assuming  a  desired  reading  on  the  order  of  1  volt  for  test 
runs,  a  drop  in  level  of  no  more  than  15.6  dB  is  indicated.  At  100  kHz, 
the  measured  attenuation  of  56.6  dB  per  meter  would  indicate  measure- 
ments could  be  made  at  least  27  cm  from  the  wedge  boundary.  This  satis- 
fies the  limit  imposed  on  r  in  Section  IV-C-1 . 
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TABLE  IV-3.  Attenuation  Results 


f(kHz) 

a(Ne/m) 

a(dB/m) 

Correlation 
Coefficient 

75.10 

7.26 

63.2 

-0.90 

100.00 

6.50 

66.6 

-0.99 

100.04 

7.39 

64.3 

-0.76 

130.00 

9.79 

85.2 

-0.92 
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V.  ERROR  ANALYSIS 

It  is  important  to  know  if  the  physical  properties  were  determined 
with  sufficient  accuracy  to  verify  Bradshaw's  model.  Bradshaw  developed 
an  analytical  equation  for  the  depression  angle  of  the  beam  formed  in 
the  substrate.  As  this  analytical  equation  was  based  on  model  results, 
and  since  the  equation  is  dependent  on  all  measured  physical  properties, 
the  equation  should  provide  a  realistic  basis  for  error  analysis.  The 
analytical  equation  is 

60  =  17.22  3°-329  9C°-772  (ow/Ps)-°-2M         («) 

where  6n  is  the  beam  depression  angle  in  degrees, 
3  is  the  wedge  angle  in  degrees, 
9  is  the  critical  angle  in  radians, 
and    p  /p  is  the  ratio  of  the  density  of  water  to  the  density  of 

saturated  sand. 
The  critical  angle  is  defined  by 

9C  =  arccos  (cw/c$)  (16) 

where  c  and  c  are  the  sound  speeds  in  water  and  sand,  respectively. 
w     s 

Using  standard  error  analysis,  and  noting  pw/ps  =  l/p  , 

deD(8,ec,Ps)=f°dB  +  !f=dec  +  f°dp      (17) 

.  5.66  9C  °-772   (l/p,)   "  °-2543°-671dB 
♦  13.2960-329   (1/PS)   "  °-25V228d6c 
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Similarly 


♦  4.37B0-329  S°-772p  0'745dp 
c    s     s 


d9c  *  ^  dcw  +  _£  dcs 
0Cw      5cs 


(18) 


w  dc  - 

s7   s 


1  del 
-   wj 


1  *  (S/Cs>' 


•1/2 


where  cw  =  1469±15  m/s  from  Table  2.  An  error  of  3  percent,  in  keeping 

with  the  worst  case  predicted  accuracy  from  Reference  10,  was  assigned  to 

sound  speed  in  saturated  sand.  Therefore,  c  =  1613£48  m/s  and  de  =  0.04 

s  c 

radians. 

Assuming  the  wedge  angle  can  be  measured  to  within  0.5  degrees,  the 
following  values  were  used  to  evaluate  equation  17. 
Ps  =  1.98  g/cm3 

3=  3° 


dp  =  0.03  g/crrT 


Therefore 


9c  =  0.4257  rad 


d9Q  =  3.64  +  0.75  +  0.16 


dS  =  0.5 


d9„  =  0.04  rad 
c 


=  4.55 
3«  for  the  above  parameters  is  15.21  degrees.  The  deviation  of  0.75 
attributable  to  sound  speed  measurements  is  roughly  5  percent,  but  is 
dwarfed  by  the  3.64  (24  percent)  error  attributable  to  the  wedge  angle 
measurement.  A  method  of  accurately  measuring  and  controlling  the 
wedge  angle  will  be  a  necessity  for  successful  model  test  runs. 
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VI.  CONCLUSIONS  AND  RECOMMENDATIONS 

A  laboratory  sized  experiment  to  test  the  predictions  of  Bradshaw's 
model  is  possible,  using  off-the-shelf  equipment  in  a  water-over-sand 
system.  Measurements  can  be  made  of  all  necessary  physical  parameters, 
and  attenuation  is  not  so  severe  as  to  preclude  measurements  at  suffi- 
cient distances  to  satisfy  model  constraints. 

Densities  and  sound  speeds  were  obtained  that  agreed  with  literature 

or  theoretical  values.  The  densities  of  water  and  saturated  sand  were 

3  3 

1.00  g/cm  and  1.98  g/cm  ,  respectively.  The  sound  speed  in  water  was 

1469  m/s;  the  sound  speed  in  saturated  sand  was  1613  m/s. 

The  use  of  a  rigid  array  of  receivers  is  highly  recommended  for  all 
measurements  within  the  sand  in  order  to  control  source-receiver  spacing. 

Model  results  are  strongly  dependent  on  the  wedge  angle.  The  sand 
will  maintain  a  flat,  sloped  surface,  but  great  care  in  the  formation  and 
measurement  of  the  wedge  angle  will  be  required. 

It  was  found  necessary  to  maintain  a  high  bleach  concentration  in 
order  to  control  gas  producing  biologic  growth.  Every  change  of  water 
should  be  accompanied  by  a  renewal  of  the  bleach  level. 
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APPENDIX  A 
SOUND  SPEED  IN  WATER  -  DATA 

TABLE  1.  Two  Receivers  Direct  Path  Data 

Source  Transducer  =  Homemade  Mylar 


Separation 
(cm) 

Time 
(10~5  sec) 

Sound  Speed 
(meters/second) 

6.8 

4.50 

1511 

11.7 

7.95 

1472 

16.7 

11.70 

1427 

19.0 

13.30 

1428 

21.4 

14.90 

1436 

26.3 

18.10 

1453 

c 

=  1454  meters/; 

second 

a 

=  32.6  (2.2%) 
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TABLE  2.  Reflected  Path,  Varied  Distance  -  Data 
Source  Transducer  =  Homemade  Mylar 


Separation 
(cm) 

Time 
(10"5  sec) 

28.8 

19.55 

38.8 

26.32 

48.8 

33.25 

58.8 

40.15 

68.8 

47.85 

78.8 

53.70 

88.8 

59.60 

98.8 

67.00 

Sound  Speed 
(meters/second) 

1473 

1474 

1468 

1464 

1438 

1467 

1490 

1475 


c  =  1469  meters/second 
a  =  14.7  (1.0%) 
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TABLE  3.  Reflected  Path,  Fixed  Distance  -  Data 
Source  Transducer  =  LC10  #2338 


Separation 

Time 

Sound  Speed 

(cm) 

(10~5  sec) 

(meters/second) 

22.0 

14.92 

1474 

22.0 

15.32 

1436 

22.0 

14.88 

1478 

22.0 

14.92 

1474 

22.0 

14.92 

1474 

22.0 

14.92 

1474 

22.0 

15.32 

1436 

22.0 

15.32 

1436 

22.0 

15.08 

c  =  1460  meters/second 
a  =  18.8  (1.3%) 

1459 
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APPENDIX  B 
SOUND  SPEED  AND  ATTENUATION  IN  SAND  -  DATA 

TABLE  1.  Sound  Speed  and  Attenuation  in  Sand 


Source  Transducer  =  Homemade  Mylar 

Frequency  =  100.04  kHz 

Source  placed  10  cm  below  sand  surface, 


Distance 

Time 

(cm) 

(10-5  sec) 

15.15 

9.04 

16.15 

10.10 

19.15 

12.25 

19.60 

12.25 

22.05 

13.80 

23.75 

15.40 

24.10 

17.55 

25.55 

15.95 

27.65 

17.00 

28.15 

19.15 

29.10 

18.05 

29.70 

18.10 

31.90 

20.20 

31.90 

21.25 

34.20 

22.85 

38.00 

no  signal 

*changed 

scale 

'*omitted 
(outside 

point 
2a) 

Amplitude  Sound  Speed 

(mV)  (meters/second) 

11.0  1676 

19.0  (18.0)*        1599 

7.0  (  7.2)         1563 

9.0  1600 

7.0  (  7.2)         1598 

3.0  (  2.8)         1542 

7.5  1373** 

4.0  (  4.0)         1602 

3.5  (  3.2)         1626 

3.0  (  2.4)         1470 

3.0  1612 

6.0  1641 

1.0  1579 

1.5  1501 

1.5  1497 


c  =  1579  meters/second 
a   =  58.6  (3.7%) 
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TABLE  2.  Sound  Speed  and  Attenuation  in  Sand 


Source  Transducer  =  Homemade  Mylar 

Frequency  =  100.00  kHz 

Source  placed  14.5  cm  below  surface  of  sand. 


Distance 

Time 

Ampl itude 

Sound  Speed 

(cm) 

(10"5  sec) 

(mV) 

(meters/second) 

10.30 

5.38 

27.0 

1614 

18.39 

11.60 

10.0 

1585 

19.15 

12.12 

42.0* 

1580 

23.71 

14.26 

4.5 

1663 

29.22 

18.94 

2.8 

1543 

*point  omitted 
for  attenuation 

c  =  1597 
a  =  44.7 

meters/ 
(2.8%) 

second 
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TABLE  3.  Sound  Speed  and  Attenuation  in  Sand 


Source  Transducer  =  Homemade  Ceramic 

Frequency  =  75.1  kHz 

Source  placed  13.0  cm  below  sand  surface 


Distance 
(cm) 

Time 
(10~5  sec) 

Ampl itude 
(mV) 

Sound  Speed 
(meters/second) 

12.2 

9.05 

25.0 

1348 

12.5 

8.00 

20.0 

1562 

23.5 

16.50 

8.0 

1424 

20.6 

14.35 

3.0 

1436 

27.4 

19.70 

2.5 

1391 

32.3 

22.35 

2.5 

1445 

c  =  1434  meters/second 
a  =  71.9  (5.0%) 
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TABLE  4.  Sound  Speed  and  Attenuation  in  Sand 


Source  Transducer  =  Homemade  Mylar 

Frequency  =  130.00  kHz 

Source  placed  7.5  cm  below  sand  surface. 


Distance 

Time 

Amplitude 

Sound  Speed 

(cm) 

(10'5  sec) 

(mV) 

(meters/second) 

12.0 

7.7 

10.0 

1558 

13.0 

8.0 



1625 

18.7 

12.0 

3.0 

1565 

23.0 

14.9 

1.8 

1544 

26.1 

16.5 

0.8 

1582 

27.3 

18.1 

0.6 

1508 

31.5 

20.2 

0.6 

1559 

32.2 

20.2 

0.6 

1594 

32.6 

20.8 

0.4 

1571 

c  = 

1567  meters/second 

a  = 

32.6 

(2.1%) 
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TABLE  5.  Sound  Speed  Gradient  Data 

Source  Transducer  =  LC10  #2338 

Frequency  =  120.05  kHz 

Source  and  receiver  fixed  at  30.0  cm  separation. 


Water 

Sand 

Time 

Amplitude 

Sound  Speed 

(cm) 

(cm) 

(10"5  sec) 

(V) 

(meters/second) 

11.0 



20.2 

5.4 

1485 

27.0 



20.2 

10.4 

1485 

30.0 

0.0 

18.5 

0.5 

1613 

29.0 

1.0 

18.6 

1.0 

1613 

27.0 

3.0 

18.6 

2.5 

1613 

25.5 

4.5 

18.6 

1.35 

1613 

24.0 

5.0 

18.6 

1.35 

1613 

23.0 

7.0 

18.6 

1.35 

1613 

21.0 

9.0 

18.6 

1.35 

1613 

20.0 

10.0 

18.6 

1.35 

1613 

17.5 

12.5 

18.6 

1.35 

1613 

15.0 

15.0 

18.6 

1.35 

1613 

12.5 

17.5 

18.6 

1.35 

1613 

10.0 

20.0 

18.6 

1.35 

1613 
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APPENDIX  C 
REFLECTION  AND  TRANSMISSION  COEFFICIENTS  -  DATA 

TABLE  1.  Surface  Reflection  Coefficient 


Direct 

Path 

Reflected 

Path 

Distance 
(cm) 

Amp. 
(V) 

Distance 
(cm) 

Amp. 
(V) 

22.5 

5.8 

48.5 

2.9 

30.0 

5.4 

37.2 

4.0 

20.0 

6.2 

40.0 

3.3* 

20.0 

6.4 

40.0 

3.0** 

R  = 

0.998 

a  = 

0.06 

*  Directivity  (.241/. 254) 
**  Directivity  (.254/. 241) 

*  and  **  measured  at  normal  incidence 


Reflection  Coeff. 

1.077 
0.9185 
1.010 
0.9880 
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TABLE  2.     Bottom  Reflection  Coefficient,   Normal    Incidence 

Distance  Directivity  Index 

Direct  =  20.0  cm  Direct  =  0.254 

Reflected  =  40.0  cm  Reflected  =  0.241 


Ampl 

itude 

(V) 

Reflection 
Coefficient 

Di  rect 

Reflected 

5.4 

0.8 

0.31 

5.7 

0.8 

0.30 

6.5 

0.8 

0.26 

7.2 

1.0 

0.29 

7.1 

1.0 

0.30 

7.1 

0.9 

0.27 

7.0 

1.0 

0.30 

7.0 

0.9 

0.27 

7.2 

1.2 

0.35 

10.0 

1.2 

0.25 

8.6 

1.1 

0.27 

R  = 

0 

29 

a   = 

0 

028 
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TABLE  3.     Bottom  Reflection  Coefficient,  Normal    Incidence 

Distance  Directivity  Index 

Direct  =  19.2  cm  Direct  =  0.254 

Reflected  =  40.4  cm  Reflected  =  0.241 


Amplitude  (V) 
Direct     Reflected 

7.2  1.1 

7.0  1.2 

7.5  1.25 

7.0  1.2 

7.1  1.0 

7.6  1.0 
7.5  1.0 
7.5  0.8 

6.5  0.8 

5.6  0.8 
6.6  0.9 
7.4  0.8 

7.4  0.95 

7.5  1.0 

7.0  1.0 

7.0  1.1 

7.5  1.2 

7.0  1.1 

7.0  0.9 
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Reflection 
Coefficient 

0.34 

0.38 

0.37 

0.38 

0.31 

0.29 

R  = 

0.31 

0.29 

a  = 

0.044 

0.23 

0.27 

0.32 

0.30 

0.24 

0.28 

0.30 

0.32 

0.35 

0.35 

0.35 

0.29 

TABLE  4.  Transmission  Coefficient,  Normal  Incidence 

Direct  Path  =  20.0  cm  in  water 

Transmitted  =  10.0  cm  in  water,  10.0  cm  in  sand 


Frequency 

(kHz) 

125 

130 

125 

120 

115 

110 

120 

125 

130 

125 

120 

115 

115 

120 

125 

130 


i  rect 

Ampl 
(vo 

itude 
Its) 
Transmitted 

Transmission 
Coefficient 

8.7 

6.8 

0.78 

9.0 

7.0 

0.78 

8.7 

5.0 

0.57 

7.8 

5.0 

0.68 

7.4 

3.6 

0.49 

6.9 

3.6 

0.52 

7.8 

5.4 

0.69 

8.7 

5.2 

0.60 

9.0 

7.6 

0.84 

8.7 

5.2 

0.60 

7.8 

5.2 

0.67 

7.4 

3.8 

0.51 

7.4 

3.8 

0.51 

7.8 

5.2 

0.67 

8.7 

5.3 

0.61 

9.0 

T  =  0 

a  =  0 

7.6 

.65 
.116 

0.84 
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